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We test approximate approaches to solving a neutrino transport problem that presents itself 
in the analysis of some accretion-disk models. Approximation # 1 consists of replacing the full, 
angular- dependent, distribution function by a two-stream simulation, where the streams are re- 
spectively outwardly and inwardly directed, with angles cos# = ±l/\/3 to the vertical. In this 
approximation the full energy dependence of the distribution function is retained, as are the energy 
and temperature dependences of the scattering rates. Approximation # 2, used in recent works on 
the subject, replaces the distribution function by an intensity function and the scattering rates by 
temperature-energy-averaged quantities. We compare the approximations to the results of solving 
the full Boltzmann equation. Under some interesting conditions, approximation # 1 passes the test; 
approximation # 2 does not. We utilize the results of our analysis to construct a toy model of a 
disc at a temperature and density such that relativistic particles are more abundant than nucleons, 
and dominate both the opacity and pressure. The nucleons will still provide most of the energy 
density. In the toy model we take the rate of heat generation (which drives the radiative transfer 
problem) to be proportional to the nucleon density. The model allows the simultaneous solution of 
Oh' the neutrino transport and hydrostatic equilibrium problems in a disk in which the nucleon density 

decreases approximately linearly as one moves from the median plane of the disk upwards, reaching 
m : zero on the upper boundary. 

PACS numbers: 4.40.Dg, 98.62.Mw 
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Hyperaccreting black holes, which have been suggested as an energy source for powering gamma ray bursts, may 
provide conditions in which the close-in part of an accretion disk reaches a region of density and temperature such 
that neutrino radiation enters the dynamics in an important way In some extreme limits, the neutrino opacity 

of the disk may become important, necessitating a solution to a neutrino transport problem. A recent paper [5| 
specifically addresses the issue of the effects of neutrino opacity in the innermost portion of such an accretion disk. 
Despite the transitory nature of the phenomena.the various time scales are such that steady flow solutions of the 
neutrino transport problem are of interest. Ref. [5j uses a simplified treatment of neutrino transport in this domain, 
a treatment that suppresses the energy dependence of the cross-sections in favor of a single parameter characterizing 
the thermal averaged opacity at an average temperature of the disk, and which further eliminates the angular degree 
^ . of freedom in the Boltzmann equation in favor of a two-stream approximation of the the angle dependence of the 
J^j ' neutrino distribution. 

The present note investigates these two approximations in domains characteristic of the accretion disks of ref.p|, 
by solving the Boltzmann equation directly, where we retain a neutrino energy spectrum and full angular dependence 
of the distribution function. The two-stream approximation to this equation, where we retain the energy dependence 
and the local temperature dependence of the opacities, is found to work quite adequately in most of the domains that 
we explore. However the further simplification of ref. |5j , in which the equation for a distribution function is replaced 
by an equation for intensities, the details of which ar7 explicated in ref. @, is not adequate in some domains that 
are of possible interest. We do not apply our results to detailed physical system, but we do sketch a technique that 
seems useful in solving the kind of neutrino transport problem that arises in these systems. We also use the technique 
to solve a toy problem of a disk in which the relativistic particles dominate the neutrino opacities and the pressure, 
but in which nucleons still provide the energy density, and in which it is assumed that the local heat generation is 
proportional to the nucleon density. Hydrostatic equilibrium can be established in this case if the nucleon density 
decreases roughly linearly with height (i.e., the coordinate transverse to the disk) and vanishes on the disk's surface. 

We consider a disk that is infinite and homogeneous in the x and y directions, with z taken as the normal coordinate. 
In appropriately scaled coordinates the surfaces of the disk are at z = ±1. The density and temperature will in general 
depend on z, with p{— z) — p(z) and T{— z) — T(z). The Boltzmann equation for the neutrino distribution function, 
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/(p, z), is now to be solved in the region < z < 1, with boundary conditions /(p, 0) = /(— p, 0), and with /(p, 1) = 

for p • z < 0. There will be a prescribed amount of heat generated per unit volume at each height in the disk, Q (z), 
and therefore a prescribed profile for the energy flux upward. The temperature function must be thus be adjusted so 
that the solution to the Boltzmann equation produces the input flux profile. 

We include both absorption and isotropic elastic scattering in the collision terms, with reaction rates, A a . s . The 
methods of solution could be extended easily enough to include the simple angular dependences of standard model 
scattering rates, but since we are trying to answer generic questions that do not depend on neutrino flavor or the 
neutron-proton mixture in the medium, we shall not include them. Experience with photon radiative transfer has 
shown that the angular dependence of Compton scattering matters little (as long at the correct weighting is included 
in the definition of the angular average of the scattering amplitude) in atmospheres that are (mathematically) similar 
to the ones that we consider p|. Likewise, as long as there is significant absorption in the opacities, inelasticity in 
scattering should not be especially consequential. 

Below we exhibit first the complete Boltzmann equation for this system, followed by two successive simplified forms. 



A. Boltzmann equation 



We define a variable x — cos(#) for the up-moving neutrinos, with cos(6>) > and x = — cos(0) for the down-moving 
neutrinos. We introduce separate notations for the distribution function in the two hemispheres in momentum space, 
f+(z,p, x) for the up- moving hemisphere and f_(z,p, x) for the down-moving hemisphere, where p is the energy (or 
momentum) of the neutrino, and where we consider absorption and isotropic elastic scattering in the collision terms. 
Then we have, 



x df+(z,p, x) = \ a [f eq (T(z),p) - f + {z,p,x)} 



dz 

f+(z,p, x) -+ 



dx' [f-(z,p,x') + f + (z,p,x')} 



(1) 



for the up-moving stream, and 



df-(z,p,x) 
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f-(z,P,x) 4 
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-Mfeq(T(z),pj - f-(z,p,x')} - 

dx' \f-{z,p,x') + f+(z,p,x')] 



(2) 



for the down-moving stream. Here f eq denotes the Fermi function with chemical potential equal to zero, and a 
temperature that depends on z. The rate constants, A a and A s are functions of p and depend on z as well, the 
latter through variation of the density of scatterers or through temperature dependence of cross-sections. Of course, 
the temperature T(z) must be adjusted so that the resulting distribution function produces the flux profile that is 
assumed, or calculated from other physics. 1 



It would be of some interest to consider cases in which the "driving chemical potential" , fl = /i e + fi p — fi n , for electron neutrinos is 
not small compared to the temperature; indeed in the supernova atmosphere, where the temperatures are a little lower than in the 
accretion disk, at the same value of the density, it is essential to do so. In this case the Fermi function f eq in and would need to 
be evaluated with chemical potential (i for neutrinos and —£l for anti-neutrinos. In a steady state calculation, then, /t(z) and T(z) (and 
hence the electron-proton ratio), must be self-consistently determined in the solution of the transport problem, simultaneously fitting 
the requirements a of given lepton number flux (zero for example) and the externally prescribed flux profile. The only treatment in the 
literature that really solves this problem in the steady flow idealization is the work of Schinder and Shapiro for the case of a hot 
newly- formed neutron star atmosphere. This work incorporates all important elastic scattering, inelastic scattering and absorption 
processes, as the present work does not. Its shortcomings are that it gives results for a rather specialized set of external parameters, 
most notably the neutrino fluxes arising from the core physics (which do not have an analogues in the disk problem) . and that it provides 
little qualitative insight into which features of the microphysics drive which features of the results. 
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B. First Simplification 

We replace the continuum angular distribution by two streams at angle cos6* = ±1/V3 to the vertical. We then 
have, 



df+(z,p) 
dz 

and, 

df-(z,p) 



X a [f eq (T(z),p) - f+(z,p)] + X s [f-(z,p) - f+(z,p)]/2 , (3) 



{) _ -X a [f eq (T(z),p)-f-(z,p)]+X s [f-(z,p)-f + (z,p)}/2, (4) 

where, X a s — \/3Aa,sj & n d where the local energy flux and energy densities are obtained from 

Flux(z) = -^-j= J p 3 dp(f+(z,p) - f-(z,p)) , (5) 

and 

PM = 4^3 / P 3 d P (f+(z,p) + f-(z,p)) • (6) 
The inside boundary condition is /_(0,p) = f+(0,p) and the outside /_(l,p)=0. 

C. Second Simplification 

If in © and the opacities were energy-independent, then we could multiply by p 3 and integrate, to get equations 
for intensity. 



dl+(z) 



X a [I eq (T{z)) - I+(z)} + ~X s [I-(z) - !+(*)] , (7) 



dz 
and 

= -X a [I eq (T{z)) + \x s [I.{z) I + (z)\ , (8) 

where I e q(T) — J dpp 3, f eq (T,p), with the inside boundary condition, 7_(0) = i+(0), and outside condition = 0. 

These are the equations discussed in ref. |6( for the case of photon radiative transfer, a solution to which, adapted to 
the neutrino case, is used in ref. Of course, for the case of neutrino transport the rate functions are dependent 
on neutrino energy as well as the temperature of the medium. Thus the set (J7J) and J5J, where we are use some sort 
of thermal averaged rate functions, represents a further, and possibly drastic, approximation to the system. 

II. INTEGRAL EQUATIONS 

We define an (absorptive) optical thickness function 



C(z,p)= / X a {y,p)dy (9) 
Jo 

The Boltzmann equations for the up-moving and down-moving distributions can now be written in integral forms, 
f + (z,p,x) = e-^^{f + (0,p,x)+ f dyx- x e^l x {x a {y,p)f eq {y,p,x) 



~ £ dx'[f + (y,p,x') + *')]])} 



+X s (y,p) -f + (y,p,x) + - dx'[f + (y,p,x') + f^(y,p,x')] , (10) 
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and 



f-{z,p,x) = e «**>/*{ J dyx- x e-^'*(\ a {y,p)U q {y,p,x) 
-K{v,p) - f-(y,P,x) + ^ J dx'[f + (y,p,x') + f-(y,p,x')]^}, 



(11) 



with f + (0,p, x) = f_(0,p,x) supplying the boundary condition in the interior. The external boundary condition is 
already incorporated into the above equations. The corresponding integral equations for the "first simplification" 
approximation are, 



f + (z,p) = e -^)(/ + (0,p) + f dye«y>ri[\ a (y,p)f eq (y,p) 



+-\ s (y,p)l-f+(y,p) + f-(y,p)} 



(12) 



and 



/_(z,p) = e«**>( [ dye~^ \\ a (y,p)f eq (y, p) 



\\ s {y,p)[~f+{y,p) + f-{y,p)]]), (13) 



where C,{z,p) = Q X a (y,p) dy. 

Beginning with an initial guess of the temperature profile and calculating the source term f eq (z,p) , we perform 
the integral in the equation for /_ in (fTTfl . taking the scattering term to be zero. We then use /+(0,p, x) = /_(0,p, x) 
as the initial condition to calculate /+ from (|10|l . This zero-order result for the distribution is used to calculate the 
scattering terms in Q and (0), the result to be used in the next iteration. In this way, if the iteration process is 
convergent, we arrive at a solution corresponding to a particular temperature distribution. Convergence problems, 
when scattering is large compared to absorption, can be avoided by feeding back, in each iteration, a scattering term 
that is, say, 75% based on the distribution from the iteration one previous to the previous one. 

The other part of our task is to adjust the temperature distribution so that the flux has the z dependence fixed (for 
now) externally. We can do this as a part of the same iteration procedure. After each iteration we calculate the flux 
from the distribution function which we then subtract (at each point z) from the desired flux function. We take these 
numbers multiplied by an appropriate constant and subtract them from the (numerically calculated) (d/dz)[T(z)]; this 
difference is taken as the new (d/dz)[T(z)], from which we construct the new function T(z) by integration, keeping 
the central temperature fixed. 

In some cases there are minor manual interventions necessary to converge to a result. In particular, for the domain 
in which the scattering is much larger than the absorption, it is necessary to ramp up the strength of the scattering 
while iterating, at a rate that keeps the errors in bounds at each stage, but that still gets the job done in a tolerable 
number of iterations. For the opacity functions that we have taken below, and for the "first simplification" form of the 
transport equation, it is surprisingly easy to get up to the range in which the scattering is six times the absorption, 
beginning in the zero'th approximation with a monotonically decreasing temperature function. For the Boltzmann 
equation itself, with its full angular dependence, this whole approach becomes much less efficient, owing to the x~ l 
factors in the exponentials, and we do not particularly recommend it as a technique, although we have used it for the 
purposes of making some of the comparisons to be summarized below. 



III. NUMERICAL COMPARISONS 



For the neutrino rate functions we have chosen two forms, which have the respective generic temperature, momen- 
tum, and baryon density dependence of the following cases: 

(I) Opacity dominated by scattering and absorption on nucleons (under conditions of sufficiently high temperature 
so that electrons are non-degenerate and the number of neutrons is approximately equal to the number of protons.) 
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X a (z,p) = f3 a n N (z)p 2 , 
X s (z,p) = (3 s n N (z)p 2 



(14) 



where tin is the nucleon number density. 
(II) Relativistic e+, e~, v, u atmosphere, 



(15) 



First we examined a solution of the "two stream" approximation in which the temperature and momentum depen- 
dences are taken into account, that is to say, the "first simplification" © and Q. We then compared with the results 
of solution of the full Boltzmann equations Q and For purposes of making the comparisons we took the flux 
to increase linearly with z. Fig. 1 compares the temperature profiles for a case with fixed temperature at the outer 
boundary, identical flux profiles, the "nucleonic" reaction rates given by (|14|) . with equal scattering and absorption, 
and common optical thickness (of the one-half disk) of 9.0 2 In this comparison and in all that follow the ordinates are 
plotted in arbitrary units; the ratios of results, once the total optical thicknesses of the disk are set, are independent 
of the absolute values of the input temperatures and densities. 
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FIG. 1: Comparison of the temperature profiles in the first two-stream simplification with that obtained from solution of the 
full Boltzmann equation, for the case of the purely nucleonic opacity function, I14H . with scattering equal to absorption. The 
flux function, surface temperature, and optical thickness equal to 9.0 are taken to be identical in the two cases. The flux is 
taken to be proportional to height, z. The units of temperature, T, are arbitrary, but the same for the two cases. The solid 
curve comes from full Boltzmann equation and the dotted curve is from the two stream approximation. 

Fig. 2 shows the same comparison for the case of purely relativistic reaction rates given in (|15[) : the agreement is 
not quite as good as that shown in fig. 1 , but together the plots show a rather impressive agreement and provide a 
strong argument for using the first form of the two-stream simplification. 

Next we address the comparison of the "first simplification" , or the two stream approximation with full temperature 
and energy dependence© and Q, and the "second simplification" , J7J and (|SJ, which subsumes all of this dependence 
in some average opacities. Popham and Narayan give an analytic solution for the latter equation for the case of z- 
independent opacities and flux that increases linearly with z, as in the two examples exhibited in figs. 1 and 2, 



(I++I-)/2 



l + f/2-(z 2 f)/2 
1 + f/2 + 1/f, 



/e(0). 



2 We define optical thickness of the half-layer in terms of the opacity functions 1141 or 1151 and the temperature profiles of the solutions 
to J1J themselves, determining a temperature-dependent mean opacity at each point, in the standard fashion, by taking the inverse 
of the appropriate thermal average of A J 1 + A^ 1 , then integrating the mean opacity (evaluated at the local temperature) to obtain the 
optical thickness. 
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FIG. 2: The same comparison as in fig.l, now for the case of purely relativistic opacity, and optical thickness of 11.0. 



where the z dependence of the source term J e , given by, 



le(0) 



(16) 



l + f/2 + l/f a -(z 2 )(f/2) 

1+ f/2 + l/fa 



le(0) 



(17) 



implicitly defines the temperature distribution. Here f and f a are the respective total and absorption depths of the 
half-layer (times the factor a/3 that enters the definition of the A functions, above). In fig. 3 we show a comparison 
of this solution with that of the "first simplified form" for the case of pure nucleonic opacity, with optical thickness 
of 11.0 and the same flux profile, this time displaying the neutrino energy densities w{z) as a function of height, in 
the two cases. The "second simplification" approximation is seen to work quite well, by the measure of the "first 
simplification" approximation, and therefore, in view of the previous comparisons, as a simulation of the full transport 
problem. 
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FIG. 3: The solid curve shows the neutrino energy density w as a function of height, as determined by the "first simplification" 
form of the equations, for the case of nucleonic opacity with equal elastic scattering and absorption rates, and optical thickness 
of 11.0, as in fig.l. The dotted curve shows the same plot for the solution l|16|l given Popham and Narayan and specialized to 
the same parameters. 

The case of purely relativistic opacity is more problematic. In fig 4. we show neutrino energy density plots, again 
for the case of equal scattering and absorption. The "second simplification" approximation here gives a greater than 
30% overestimate of the (central energy density)/ (flux) ratio. But more significant, perhaps, are the shapes of the 
curves; with the same flux profile, the "first simplification" approximation has a quite different energy density profile 
from the "second simplification" result. In particular, note that the negative of the gradient of the energy density, in 
the outer region, is a decreasing function of z in the entire region z > .5 in our treatment that takes into account the 
energy and temperature dependences, while it is positive throughout in the case in which averaged opacities are used. 
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FIG. 4: The same as fig. 3, but for purely relativistic opacity, with optical thickness 11.0 as in fig. 1 . 



When we plot the T(z) 4 instead of the energy density of the neutrinos we find a shape that differs appreciably 
from that shown in fig. 4 only in a region extremely near the surface (in fact, a region beginning above the average 
height for the last neutrino scattering). Thus the gradient of the pressure of all of the relativistic particles, taken 
together, declines steadily as we move upwards in the upper parts of the disk. Hydrostatic equilibrium requires that 
the total pressure gradient be proportional to zp(z), where p(z) is the total energy density, since the operative tidal 
gravitational field in the disk is proportional to z. Therefore if relativistic pressure were the whole story, the density 
would have to be a fairly rapidly decreasing function of z, in contrast to the uniform density assumption that appears 
to be used in In any of the regimes that we here consider, the total energy density is dominated by nucleons. But 
both pressure and neutrino opacities are dominated by relativistic particles when the temperature is such that the 
ratio of the number of relativistic particles to the number of nucleons is large. 

One can therefore find regimes in which the relativistic particles predominate in both the pressure and opacity 
calculations while the mass density is still dominated by nucleons. We might achieve hydrostatic equilibrium in such 
a system by having a nucleon density that decreases significantly from the center to the top. For example, consider 
a nucleon density that is proportional to (1 — \z\). Now assume as well that heat production is proportional to this 
density. Then we need to solve the neutrino transport problem, taking a flux profile proportional to z — z 2 Jl (for 
z > 0). The resulting relativistic pressure profile is shown in fig. 5 for the case of (optical thickness=11.0). 3 




FIG. 5: The solid curve shows the total pressure (in arbitrary units) of relativistic particles as a function of height, as determined 
by the "first simplification" form of the equations, for the case of relativistic opacity with equal elastic scattering and absorption 
rates, where the flux is now taken as proportional to z — z 2 /2. The optical thickness is 11. The dotted curve shows the function, 
a + b(z 2 /2 — 2 3 /3) that most nearly fits the pressure curve 



3 The total relativistic pressure is the sum of the neutrino pressure and that of the e^ particles and the photons, the latter two computed 
from the local temperature; for the neutrino contribution we take the local energy density, divided by 3, from the solution of the radiative 
transfer problem. The shape of the energy density curve is the same as that of the total relativistic pressure except in the very outermost 
region. 



8 



On the same curve we have drawn the profile of the pressure required for hydrostatic equilibrium for the case in 
which the total energy density is proportional to (1 — \z\), giving P = const x (z 2 /2 — |z| 3 /3). It is clear that with 
a small redistribution of the nucleons, in this model we could simultaneously achieve our three goals: (a) Solution of 
the radiative transfer problem; (b) hydrostatic equilibrium; and (c) matching the energy flux to a local heating rate 
that is proportional to the nucleon density. 

In the case of accretion disks that are more or less uniform in density it appears that the temperature will never 
be quite high enough, at a given nucleon density, to lead to the complete dominance of pressure and opacity by the 
relativistic part of the medium |5|- But the above example does suggest the possibility of stable relativistic outlying 
regions, still well under the height of (average) last neutrino scattering, in which nucleon densities decrease markedly 
with height. The surface of the disk cannot be as simple as in this idealized model and one can easily imagine 
instabilities that lead to escape of some of the surface plasma. Thus it might be possible to achieve an electron- 
positron plasma in the region both near the surface and the axis that is more dense than that coming from neutrino- 
neutrino collisions in the vicinity of the axis, the latter a source that has been discussed in the literature 0. 

The author is indebted to Andrei Beloborodov and to Chris Fryer for interesting conversations. 
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